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We prove the theorem: The second order quasi-hnear differential operator as a second rank 
divergence free tensor in the equation of motion for gravitation could always be derived from the 
trace of the Bianchi derivative of the fourth rank tensor, which is a homogeneous polynomial in 
curvatures. The existence of such a tensor for each term in the polynomial Lagrangian is a new 
characterization of the Lovelock gravity. 
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The Bianchi differential identity involving the purely 
anti-symmetric derivative of a derivative {D^ = 0) 
was famously interpreted by John Wheeler [l[ as the 
statement of the fact, boundary of boundary is zero. 
The familiar examples of it are curl of gradient and 
divergence of curl being zero. The former signifies 
a scalar while the latter a vector field. However its 
contraction in the two cases is vacuous and hence 
docs not lead to a non-trivial statement I^'. When we 
go beyond vector to a tensor field, it becomes interesting. 

Gravity is the universal force which means it links to 
all particles unmindful of their mass being non-zero or 
zero. It is its linkage to massless particles which leads to 
the profound realization that it could only be described 
by curvature of spacetime [H, This means the dy- 
namics of gravity resides in spacetime curvature which 
must fully and entirely determine it. That is the gravita- 
tional dynamics must follow from the curvature which is 
described by the fourth rank Riemann curvature tensor 
(it is defined as A!;^. — A\i = R^mikA™, a generalized 
"curl"). It involves second and first derivatives of the 
spacetime metric, gat- The Bianchi identity is given by 
the anti-symmetric derivative of R^mik, 



(1) 



If the gravitational dynamics has to follow from the cur- 
vature, it has to follow from this identity which is the 
only available geometric relation. The only thing we can 
do to it is to contract on the available indices which docs 
lead, unlike for scalar and vector, to a non- vacuous rela- 
tion, 

(^""bia ~ 0, Gab ~ Rat ~ -^^9ah (2) 



where Rah is the Ricci tensor, the contraction of Rie- 
mann, while R is the trace of Ricci. Now the trace (con- 
traction) of the Bianchi identity yields a non-trivial dif- 
ferential identity from which wc can deduce the following 
relation 



Gah — KTab ~ ^Qab, T°' 







(3) 



where Tab is the second rank symmetric tensor with 
vanishing divergence and k and A are constants. The left 
hand side of the equation is a second order differential 
operator on the metric gab- For this equation to describe 
dynamics of gravity, the tensor Tab should describe the 
source/charge for gravity which should also be universal. 
It should be shared by all particles and hence Tab should 
represent energy momentum distribution. Thus we ob- 
tain the Einstein equation for gravitation which entirely 
follows from the spacetime curvature. We have however 
two constants of which one k is to be determined by 
experimentally measuring the strength of the force and 
is identified with Newton's constant, k — SttG/c^. 
Why is there new constant A which though arises in 
the equation as naturally as the energy momentum 
tensor, T^h? It is perhaps because of the absence of 
fixed spacetime background which exists for the rest of 
physics and the new constant may be a signature of this 
fact. It is the universal character of gravity which makes 
spacetime itself dynamic. The force free state would 
however be characterized by homogeneity and isotropy 
of space and homogeneity of time which will in general 
be described by spacetime of constant curvature and not 
necessarily of zero curvature. The new constant A is the 
measure of the constant curvature of spacetime and it 
identifies the most general spacetime for force free state. 
It may in some deep and fundamental sense be related 
to the basic structure of spacetime. 
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[1] For a vector field, vanishing of Bianchi derivative means the field 
is a curl of gauge potential;i.e. Fif^ = fe — ^fc i = —Fki- This 
trivially implies F}^^, = and consequently we can write F}^ = J' 
with J* = 0. This is indeed the equation of motion for the gauge 
field with J* being identified as the current vector. 



We also know that the complete contraction of Rie- 
mann gives the scalar curvature, i?, the Einstein-Hilbert 
Lagrangian which on variation leads to the divergence 
free Einstein tensor. Gab, and subsequently the Einstein 
equation. We thus have the two distinct but equivalent 
derivations for the gravitational dynamics. The former 
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is simply driven by the geometry while the latter is in 
the spirit of every dynamical equation following from 
an action. It is always possible to write an action con- 
structed from Riemann curvature for higher derivative 
gravity and derive the corresponding equation of motion. 
Similarly, is it possible to derive an analogue of Gab, a 
divergence free differential operator from the Bianchi 
derivative of the higher order curvature polynomial? 
This is the question we wish to address and show that 
the answer is in affirmative. It would give yet another 
characterization of the Lovelock gravity. 

We believe that gravitational dynamics to follow from 
the spacetime curvature should be a general principle 
which should be true in general for higher order theo- 
ries as well. So comes the question of going beyond the 
linear order in Riemann. Let us consider the quadratic 
tensor, 

Tlabcd = RahmnRcd™'^ + Oi:R[a"^ R,b],ncd + PRRabcd (4) 

where a, f3 are constants. We consider the Bianchi deriva- 
tive, TZablcd;e] which On contraction gives 

where TZac — g'"^TZabcd and TZ ~ g''''flab- It turns out 
that for a = 4 and /3 = 1, we obtain 

= {-H/' + ^LGB6el;c. (6) 

That is 

-^""[cd^e] ' liLGB5el,c = -H,^, = 0. (7) 

The tensor Hab is divergence free, Ha\ = 0, and is given 
by 

liab — Z(KKab — ^bm ~ ^n. Uarnbn 

+ Ra™^^Rbmnl) — -^LcBgab- (8) 

It results from the variation of the well-known Gauss- 
Bonnet Lagrangian Lq b = Rabcd ~ ^Rab + R^ where we 
have written R^bcd ^ RabcdR""^"^ ■ That is, we can write 

Hab = ^.Tlab - \ngab (9) 

where TZ ~ Lqb- 

Though Hab can be written in terms of TZabcd but it 
doesn't follow directly from it as Gab does from Rabcd- 
We note that Bianchi derivative vanishes only for the 
Riemann curvature signifying the fact it can be written 
as a generalized curl of a vector. No other tensor will 
have vanishing Bianchi derivative. However to find the 
analogue of Gab, we donot require vanishing of Bianchi 



derivative of the quadratic tensor TZabcd but instead van- 
ishing of its trace would suffice. We see that even the 
trace of Bianchi derivative does not vanish but is equal to 
\{Lgb5c'^):c- It suggests that the curvature polynomial 
should also include term involving its own trace which 
would make no contribution in the linear case. Before 
we do that, let us write TZabcd and Hab for a general or- 
der n in the Lovelock polynomial @. In the context 
of higher derivative gravity theories, a unified scheme of 
writing Lagrangian is given [3] as an invariant product, 
Q'^'^'^'^ Rabcd where Q"'^'^'^ has the same symmetry proper- 
ties as the Riemann tensor. It is constructed from metric 
and Riemann curvature and has vanishing divergence, 
Qabcd _ g p^j. -j-j^^, Lovclock Lagrangian, it is required 
to be a homogeneous function of the Riemann curvature. 
Since TZabcd is also a homogeneous polynomial in Rie- 
mann curvature for a given order, it is hence possible to 
write 

TZabcd = Qab"^^Rmncd (10) 

where Qabcd is in general given for n-th order by 0], 

Qab _ ra6ai6i...a„6„ p cidi p c„(i„ ("1 1 \ 

cd - Ocdcicii...c„d„-n^aifci ---Ra^b^ ■ (.ttj 

Let us verify this for the quadratic Gauss-Bonnet case 
where Qabcd could explicitly be written as 0, 

Qabcd = Rabcd " 2i?a[c5d]h + 2Rb[cgd]a + Rga[cgd]b- (12) 

It is easy to see that when it is substituted in Eqn. (|10p , 
we obtain the Gauss-Bonnet TZabcd as given in Eqn. ^ 
with a ~ i, (3 ~ 1. We could thus write TZabcd by using 
Eqs pT|) and (fTU]) for any term in the Lovelock polyno- 
mial and the corresponding Hab for the 72-th polynomial 
will be given by 

Hab = nTZab - ^TZgab (13) 

which will for the linear case n = 1 be the Einstein ten- 
sor. Gab- In this case, the trace of the Bianchi derivative 
vanishes because the Riemann tensor satisfies the Bianchi 
identity. However Hab as defined above is the analogue of 
the Einstein tensor for n-th order polynomial and hence 
the generalized Einstein tensor is divergence free. As 
trace of Rabcd is the Einstein-Hilbert Lagrangian simi- 
larly trace of TZabcd gives the Gauss-Bonnet and higher 
order Lovelock Lagrangian. Note that Eqn. ([S]) holds 
good in general for any order n of the polynomial TZabcd- 
We can thus write in general for Eqn. 

TZ^\cd:e] - "^{TZS/').., = --H,% = (14) 
n n ' 

where TZ is the corresponding n-th order Lagrangian. For 
71 = 1, the second term on the left vanishes indicating 
the Bianchi identity and Hab = Gab- For n > 1, the 
trace of Bianchi derivative does not vanish unless the 
trace TZ is also included as indicated by the second term 
which then leads to the required divergence free tensor. 
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Hab- Apart from the quadratic Gauss-Bonnet case, we 
have also verified the above relation for 71 = 3 cubic 
polynomial. It is interesting that the order of polynomial 
appears in the definition of the corresponding curvature 
polynomial tensor yielding the second rank symmetric 
divergence free tensor. 

We now turn to redefining the curvature polynomial 
which also includes its trace such that trace of the Bianchi 
derivative vanishes directly giving the required differen- 
tial operator in terms of divergence free Hab- We thus 
write 



n.brd T^ahrd 



1 



'R-igacgbd-gad.gbc) (15) 



where n being the order of the polynomial and d the 
dimension of spacetime under consideration. Then 



■cd 







(16) 



which verifies Eqn. (|13[). In terms of J- abed, Hab is given 
by 



niJ^ab - gab) = H, 



ah 



(17) 



where Fab = g J^acbd and J" = g Tab- It is interesting 
to note that Tab is the analogue of the Ricci, Rab and 
Hab of the Einstein, Gab with the order n of the polyno- 
mial TZabcd which is the analogue of the Riemann, Rabcd- 
Thus Hab is obtained from J- abed in the same manner as 

Gab is from Rabcd, and for n = 1, Tabed = T^abed = Rabed- 

AU this however happens only for the specific Lovelock 
coefficients in the polynomial and thereby determining 
the Lovelock polynomial. 

Note that Hab defines the general conserved (Einstein) 
tensor for any order n of the homogeneous curvature 
polynomial and let us take its trace which would give 



d-2n 
n{d-2) 



n. 



(18) 



For the usual Einstein gravity, n = 1, it is = 7?. 
with d > 2 and the Einstein tensor vanishes for c? = 2. 
Thus we have the general result that d > 2n for n-th 
order polynomial and Hab vanishes for d = 2n. It is 
the variation of TZ (and so also T) gives Hab which 
vanishes for d = 2n but not TZ. If we take J- which is 
proportional to the Lovelock Lagrangian TZ, then the 
Lagrangian itself vanishes for d = 2n. It is remarkable 
that the existence of an analogue of the Einstein tensor 
in general - a conserved tensor which is non-zero only 
for d > 2n, uniquely determines the Lovelock polynomial. 

By redefining the curvature polynomial we have been 
able to derive the divergence free differential operator 
for the gravitational dynamics. This is what we had set 



out to do and it could be stated as follows: 

Theorem: The second order quasi-linear differential 
operator as a second rank divergence free tensor in 
the equation of motion for gravitation could always 
he derived from the trace of the Bianchi derivative of 
the fourth rank tensor, Tabcd, which is a homogeneous 
polynomial in curvatures. The trace of the curvature 
polynomial is proportional to the corresponding term 
in the Lovelock action and corresponding to each term 
in the Lovelock Lagrangian, there exists a fourth rank 
tensor which is a new characterization of the Lovelock 
Lagrangian. 

It is the requirement of quasi-linearity (linear in sec- 
ond derivative) of the equation of motion which singled 
out the Lovelock polynomial with specific coefficients. In 
our case it is replaced by the requirement that the fourth 
rank tensor which is homogeneous in curvatures yields 
a divergence free second rank tensor through the trace 
of its Bianchi derivative. This automatically ensures the 
quasi-linearity of the equation. For non Lovelock action, 
there would not exist a fourth rank tensor with this 
property. 

In a sense, the Riemann curvature could be looked 
upon as a Bianchi potential giving the Einstein tensor 
for the Einstein gravity while Tabcd is the Bianchi 
potential giving Hab for the Lovelock dynamics. Thus 
each term in the Lovelock polynomial has a potential 
tensor Fabed given by Eqs dTO]), (HI]) and It 
is remarkable that there exists Bianchi potential for 
each term in the Lovelock Lagrangian. Non existence 
of potential characterizes all other actions like f{R) 
gravity. The existence of Bianchi potential thus becomes 
a distinguishing feature for the Lovelock action. In this 
context it may be noted that the requirement that both 
Palitini and metric action give the same equation of 
motion also picks up the Lovelock action @. We have 
thus three distinct properties (namely quasi-linearity of 
equation of motion, equivalence of Palitini and metric 
formulation and existence of Bianchi potential) which 
characterize the Lovelock action. I believe that there 
should exist a thread knitting them and it would be 
interesting to probe that. 

However, the relevance of the order in the Lagrangian 
for the gravitational dynamics depends upon the 
spacetime dimension. For instance, for d < 5, the 
quadratic Gauss-Bonnet term makes no contribution 
to the equation of motion and similarly the cubic term 
becomes relevant only for d > 6. We would however like 
to emphasize that for d > 4, Hab must be included along 
with Gab for the classical dynamics of gravitation. In 
the ultra violet limt of the theory signifying super strong 
gravitational field, it is pertinent to include higher order 
curvature effects. If we continue to have a well defined 
evolution of the field, the equation must be quasi-linear. 



4 



This will inevitably and uniquely lead to Lovelock 
polynomial and higher dimension. Further one loop 
correction in string theory generates the Gauss-Bonnet 
term Q. That is, strong field limit of classical gravity 
and one loop quantum correction seem to share the 
same Gauss-Bonnet ground. It could therefore be 
envisioned that the classical limit to quantum gravity 
is perhaps via the Lovelock gravity and the relevant 
order in the polynomial (loop correction) being given 
by the spacetime dimension under consideration. The 
Gauss-Bonnet and higher orders may therefore represent 
a intermediary state between the classical Einstein 
gravity and quantum gravity Q ■ 

Gravity is an inherently self interactive force and the 
self interaction could only be evaluated by successive 
iterations. The Einstein gravity is self interactive but it 
contains only the first iteration through the square of 
first derivative in Ricmann curvature. The question is 
how do we stop at the first iteration? The second iter- 
ation would ask for a quadratic polynomial in Riemann 
curvature which should give the corresponding term in 
the equation of motion. Thus the quadratic tensor, TZabcd 
as given in Eqn. (|4]) with specific coefficients, will alone 
meet the requirement (quasi- linearity of the differential 
operator) for inclusion of the second iteration. However 
its effect in the equation of motion can be felt only for 
dimension, o? > 4, and hence we have to go to higher 
dimension for physical realization of the second iteration 
of self interaction 0, [E]- It is remarkable that even 
classical dynamics of gravity asks for dimension > 4. As 
two and three dimensions were not big enough for free 
propagation of gravity, similarly four dimension is not 
big enough to fully accommodate self interaction dynam- 
ics of gravity. Then the most pertinent question is where 
does this chain end? Let us envision the AdS/CFT-like 
scenario where the 3-brane ((3 -I- l)-spacetime) forms 
the boundary enclosing the higher dimensional bulk 
spacetime. If matter fields remain confined to 3-brane, 
the bulk would then be free of matter and hence it 
would be maximally symmetric (homogeneous and 
isotropic in space and homogeneous in time). It is 
then a spacetime of constant curvature, dS/AdS, with 
vanishing Weyl curvature. There is no free gravity to 
propagate any further and hence the chain stops at 
the second level at least in this particular construction. 
Whether this is a generic setting or not is however an 
open question? It may be noted here that gravitational 
equation with inclusion of both Gab and Hab for empty 
space in higher dimension admits dS / AdS as solution [J] . 



In what dimension should matter live is to be entirely 
determined by the matter dynamics. If the matter fields 
are gauge vector fields described by 2-form, the confor- 
mal invariance (universal scale change, gab ^ .f^9ab) 
of their dynamics will dictate that they live in four 
dimension. This general principle is always obeyed by 
the matter field dynamics unless a scale in terms of mass 
etc. is introduced by spontaneous symmetry breaking. 
The symmetry breaking is indicative of theory being 
incomplete and it is hoped that the complete theory 
would restore the conformal invariance. It is therefore 
reasonable to take that matter remains confined to 
3-brane. This view is also supported by the string 
theory where open strings have their endpoints on the 
brane indicating residence of matter there. 

The most fundamental question is how do we know 
there exist higher dimensions and if so why do we not 
see them? The existence of dimension can only be 
probed by a physical interaction. All our probes are 
matter field forces like electromagnetic, which remain, as 
argued above, confined to 3-brane and hence they cannot 
fathom higher dimension. Since gravity is universal 
and hence it cannot be confined entirely to the brane 
and can propagate in higher dimension. Above we have 
argued quite convincingly that there is a strong case for 
higher dimension for physical realization of the second 
iteration of self interaction. It is a purely classical 
motivation (we have also elsewhere 0, Q given a couple 
of more classical arguments) while higher dimension is 
a natural arena for string theory. The only way higher 
dimension could thus be probed is by a very high energy 
pure gravitational experiment. This is what we have not 
yet been able to do and hence the question remains open. 

We had set out to establish the general principle that 
gravitational dynamics resides in spacetime curvature 
and hence it should always and entirely be driven by 
spacetime geometry. This we have shown by deriving 
the quasi-linear differential operator for the equation of 
motion for the Lovelock gravity. In particular, we have 
found a new geometric characterization of the Lovelock 
gravity in existence of Bianchi potential for each term 
in the polynomial. This is indeed an interesting general 
property. 
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